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Abstract 

Accretion disk winds browing off perpendicular to a luminous disk are examined in the framework 
of fully special relativistic radiation hydrodynamics. The wind is assumed to be steady, vertical, and 
isothermal. Using a velocity-dependent variable Eddington factor, we can solve the rigorous equations 
of relativistic radiative hydrodynamics, and can obtain radiatively driven winds accelerated up to the 
relativistic speed. For less luminous cases, disk winds are transonic types passing through saddle type 
critical points, and the final speed of winds increases as the disk flux and/or the isothermal sound speed 
increase. For luminous cases, on the other hand, disk winds are always supersonic, since critical points 
disappear due to the characteristic nature of the disk gravitational fields. The boundary between the 
transonic and supersonic types is located at around F c ~ 0.1(e +p) / (pc 2 )/^ c , where F c is the radiative flux 
at the critical point normalized by the local Eddington luminosity, (e + p) / (pc 2 ) is the enthalpy of the gas 
divided by the rest mass energy, and j c is the Lorentz factor of the wind velocity at the critical point. In 
the transonic winds, the final speed becomes 0.4-0. 8c for typical parameters, while it can reach ~ c in the 
supersonic winds. 

Key words: accretion, accretion disks — astrophysical jets — gamma-ray burst — radiative transfer 
— relativity 
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1. Introduction 

Mass outflows - accretion disk winds - from a lumi- 
nous disk are observed in various active objects, such 
as cataclysmic variables (CVs), supersoft X-ray sources 
(SSXSs), microquasars (/xQSOs), broad absorption line 
quasars (BAL QSOs), and so on. Radiatively driven wind 
emanating from an accretion disk is a clue to the forma- 
tion mechanism of astrophysical jets and winds in these 
objects (see Kato et al. 1998, 2007 for a review of accre- 
tion disks). 

So far, radiatively driven outflows from a luminous 
disk have been extensively studied by many researchers 
(Bisnovatyi-Kogan, Blinnikov 1977; Katz 1980; Icke 1980; 
Melia, Konigl 1989; Misra, Melia 1993; Tajima, Fukue 
1996, 1998; Watarai, Fukue 1999; Hirai, Fukue 2001; 
Fukue et al. 2001; Orihara, Fukue 2003), and by numerical 
simulations (Eggum et al. 1985, 1988; Kley 1989; Okuda 
et al. 1997; Kley, Lin 1999; Okuda, Fujita 2000; Okuda 
2002; Okuda et al. 2005; Ohsuga et al. 2005; Ohsuga 
2006). In almost all of these studies, however, the lumi- 
nous disk was treated as an external radiation source (op- 
tically thin approximation), and the radiation transfer in 
the flow was not solved. Although radiation hydrodynam- 
ical equations were solved in the numerical simulations, 
the methods were rather limited; e.g., the equilibrium be- 
tween gas and radiation was assumed, the flux-limited dif- 
fusion approximation was adopted, the flow velocity was 



Present address: Center for Computational Physics, University 
of Tsukuba, Tennoudai 1-1-1, Tsukuba, Ibaraki, 305-8577 



subrelativistic on the order of 0.1c, or the optically thick 
to thin transition was not properly treated. Up to now, 
in relation to accretion disk winds no one solved the fully 
relativistic radiation hydrodynamical equations. 

Recently, radiation hydrodynamical mass outflows have 
been examined for the first time in the fully relativistic 
cases for the plane-parallel case (Fukue 2005, 2006; Fukue, 
Akizuki 2006; see also Akizuki and Fukue 2007 for the 
spherical case). In these studies, however, the gas pres- 
sure was ignored for simplicity, although the relativistic 
radiation hydrodynamical equations were solved. In this 
paper, we thus take into account the gas pressure, and ob- 
tain the relativistic radiation hydrodynamical winds brow- 
ing off from a luminous flat disk within the framework of 
a fully special relativistic regime. 

In the next section we describe the basic equations in 
the vertical direction, and examine critical points. In sec- 
tion 3 we obtain transonic winds as well as supersonic 
ones. The final section is devoted to concluding remarks. 

2. Basic Equations and Boundary Conditions 

Let us suppose a luminous flat disk, inside of which 
gravitational or nuclear energy is released via viscous heat- 
ing or other processes. The radiation energy is trans- 
ported in the vertical direction, and the disk gas, itself, 
also moves in the vertical direction due to the action of 
radiation pressure (steady plane-parallel approximation). 
We do not consider the rotation of the gas. As for the or- 
der of the flow velocity v, we consider the fully relativistic 
regime, where the terms are retained up to the second or- 



der of (v/c). As for the gravity, on the other hand, we 
adopt the pseudo-Newtonian approximation (Paczyhski, 
Wiita 1980), since we do not consider the region very 
close to the Schwarzschild radius. For simplicity, in the 
present paper, we assume that the gas is isothermal, be- 
cause we focus our attention on the general properties of 
transonic disk winds driven by disk radiation fields under 
relativistic radiation hydrodynamics. We further assume 
the radiative equilibrium and use the gray approximation. 
Finally, in order to close moment equations, we adopt the 
velocity- dependent variable Eddington factor proposed by 
Fukue (2006). 

2.1. Basic Equations 

Under these assumptions, the radiation hydrodynamic 
equations for steady vertical (z) winds are described as 
follows (Kato ct al. 1998, 2007; Fukue 2006). 

The continuity equation is 
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where p is the proper gas density, u the vertical four ve- 
locity, J the mass-loss rate per unit area, and c the speed 
of light. The four velocity u is related to the proper three 
velocity v by u — jv/c, where 7 is the Lorentz factor, 

7 = VT^u2 = l/^l_( v / c )2. 

The equation of motion is, within the present approxi- 
mation. 
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where M is the mass of the central object, R = \/r 2 + z 2 , r 
being the radius, r g (= 2GM/c 2 ) the Schwarzschild radius, 
e the gas internal energy per unit proper volume, p the gas 
pressure measured in the comoving frame, K a b s and K sca 
the absorption and scattering opacities (gray), defined in 
the comoving frame, E the radiation energy density, F 
the radiative flux, and P the radiation pressure observed 
in the incrtial frame (Fukue 2006; Kato et al. 2007). The 
first term in the square bracket on the right-hand side of 
equation (2) means the radiatively-driven force, which is 
modified to the order of it 2 , whereas the second term is the 
radiation drag force, which is also modified, but roughly 
proportional to the velocity. Compared with the previous 
researches (Fukue 2006; Fukue and Akizuki 2006), we have 
added the pressure gradient force, the second term on the 
right-hand side of equation (2), and the related factor, 
pc 2 /(e + p) to the third term. 

The radiative equilibrium condition is, in the inertial 
frame, written as 



= j - K^bscEj 2 - K ahs cPu 2 
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where j is the emissivity defined in the comoving frame 
(Fukue 2006; Kato et al. 2007). 

For radiation fields, the zeroth-moment equation be- 
comes 



7 



[F(l + 2u 2 ) 7 - (cE + cP) 7 2 u] , 



(4) 



where we use equation (3) for the second equality. The 
first-moment equation is 

dP p7 

~T = — b'( M /7) - K ahs F + K ahs cP(u/"f) 
dz c 
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where we use equation (3) for the second equality (Fukue 
2006; Kato et al. 2007). 

In order to close moment equations for radiation fields, 
we need some closure relation. Instead of the usual 
Eddington approximation, we here adopt a velocity- 
dependent variable Eddington factor /(/3), 



P Q = f(f3)E 



(0) 



in the comoving frame, where Pq and Eq are the quantities 
in the comoving frame, and (3 — v/c. If we adopt this 
form (6) as the closure relation in the comoving frame, 
the transformed closure relation in the incrtial frame is 

cP (1 + u 2 - fu 2 ) = cE (/ 7 2 - u 2 ) + 2F 1U (1 - /) , (7) 
or equivalcntly, 

cP (1 - f(P) =cE(f- (3 2 ) + 2FI3 (1 - /) . (8) 
As a form of the function /(/?), we adopt the simplest one: 



(9) 



for a plane-parallel geometry (Fukue 2006; Fukue, Akizuki 
2006; cf. Akizuki and Fukue 2007 for a spherically sym- 
metric geometry). 

Using continuity equation (1) and equation of motion 
(2), after some manipulations, we obtain the so-called 
wind equation, 

du u ( GMz 
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or equivalcntly, 
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Here, ct is the constant isothermal sound speed, defined 

by 



2 C P 
C T = 

e + p 



(const.). 



(12) 



Under the present isothermal assumption, the factor, 
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e + p 

is also a constant parameter. 

Eliminating p and E with the help of equations (1) and 
(7), equations (11), (4), and (5) become 
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For the optical depth r, defined by dr s 
we have 

^sna. 1 



j K abs 

dz C 



-(/Cabs + Ksca)P<fe) 



(17) 



Finally, using the nondimensional variables, such as z = 
z/r g , (3 = v/c, F = F/[L E /(A- ! rrl)], P = cP/[L E /{4xr*)], 
where Le [= AttcGM / (K a i, s + K sca ,)] is the Eddington lu- 
minosity, basic equations (14)-(17) are normalized as 
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18) 
(19) 
(20) 
(21) 



where c| = c^/c 2 = p/(e+p), 6 2 = pc 2 /(e+p), and J = 
C 2 J/[L E /(4^r 2 )]. 

We solve equations (18)— (21) under regularity condi- 
tions at the sonic point and appropriate boundary con- 
ditions at the moving surface with a variable Eddington 
factor (9). 

2.2. Regularity Conditions and Boundary Conditions 

As is well known, equation (18) has critical points, 
where the denominator and numerator vanish, simulta- 
neously. Hence, transonic solutions, which passes through 
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Fig. 1. Loci and types of critical points as a function of 
the radiative flux _F C at the critical point for several values of 
the isothermal sound speed fix- Filled symbols represent the 
saddle type, whereas open symbols denote the center type. 
The values of ct are attached on each symbols. The other 
parameters are fixed as 5 = 1, r c = 1, and P c = 0.1. 

critical points, satisfy the regularity conditions at the crit- 
ical point. In order to obtain transonic wind solutions, in 
this paper, we first search the location of critical points, 
examine the types of critical points, and calculate the 
transonic solutions from the critical point both inward 
and outward directions. 

In contrast to a simple spherical flow under the gravity 
of the central object, the gravitational field of the disk 
wind does not monotonically decrease, but increases at 
first and then decreases (Fukue 2002). As a result, there 
may appear multiple critical points. In figure 1, we show 
some typical examples. 

In figure 1, the height z c of critical points are plotted 
as a function of the radiative flux F c at the critical point 
for several values of the isothermal sound speed ct . Filled 
symbols represent the saddle type, through which physical 
solutions can path, whereas open symbols denote the cen- 
ter type, where no solution can path (Kato et al. 2007). 
The values of ct are attached on each symbols. The other 
parameters are fixed as S = 1, r c = 1, and P c = 0.1. 

As is easily seen in figure 1, in less- luminous cases there 
appear two positions for the same values of P c . Of these, 
the upper point is usually a saddle type, whereas the lower 
one is a center type. On the other hand, in luminous cases 
there is no critical point. Roughly speaking, two critical 
points exist when i5 2 7 c i^ c < 0.1, or 



W(4^r 2 ) 



<0.1 



e+p 1 

pc^ 7c 



(22) 



where F c is the radiative flux at the critical point, (e + 
p) I (pc 2 ) is the enthalpy of the gas divided by the rest mass 
energy, and 7 C is the Lorentz factor of the wind velocity 
at the critical point. 

In addition to the above regularity conditions, at the 
wind top of r = we must impose the boundary condi- 
tions, which is different from the usual ones for a static 
photosphere, since the wind top moves upward at a rel- 
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ativistic speed (Fukue 2005). In the case of a "moving 
photosphere", due to relativistic aberration and Dopplcr 
effect (cf. Fukue 2000), the boundary conditions imposed 
for the radiation quantities become 



CPs 

F s 



2 + 6& + 6& 2 

3 + 8/3 s + 3/3 h 2 ' 



(23) 



where the subscript 's' denotes the quantitcs at the wind 
top, and (3 a is a final speed at the wind top. 

3. Relativistic Accretion Disk Winds 

In this section we show typical solutions of relativistic 
radiation hydrodynamical accretion disk winds for tran- 
sonic and supersonic cases. 

In order to obtain a transonic solution, by numerically 
solving equations (18)-(21), we must specify the initial 
values of r, and z c , (3 C , F c , P c , r c at the critical point for 
a given set of parameters of ct, 8, and J. However, there 
are three restrictive conditions; two regularity conditions 
at the critical points, and one boundary condition at the 
wind top. Hence, we have six freedoms. In order to see 
the effect of the radiative force and the gas pressure, in 
the present paper we set r = 3, r c = 1, S = 1, and J = 0.1. 
If we give the values of F c and ct, the values of z c , (3 C , 
and P c are automatically determined. 

Several examples of transonic solutions are shown in 
figure 2. Physical quantities arc normalized in terms of 
the speed of light c, the Schwarzschild radius r g , and the 
Eddington luminosity Le [= 4ncGM/ (« a bs + K sca )} ; the 
units of F and cP are LE/(47rrg). 

In figure 2 we show the flow velocity (3 (thick solid 
curve), the radiative flux F (dashed one), the radiation 
pressure P (dotted one), and the optical depth r (solid 
curve), as a function of the height z for f = 3, 5 = 1, r c = l, 
F c = 1, J = 0.1, and ct =0.1, 0.3, 0.5. Critical points are 
marked by open circles. 

As is seen in, e.g., figure 2a, the wind velocity slightly 
decreases at first, because of the existence of the center- 
type critical point, then increases to pass through the 
saddle- type critical point, and finally reaches the final 
speed at the wind top of r = 0. The radiative flux F 
slightly decreases, as the height increases. In the relativis- 
tic radiation flow, the radiative flux does not conserve, but 
decreases, since the radiation field acts to accelerate the 
gas. As for the effect of the gas pressure, the velocity fields 
entirely increase as the isothermal sound speed increases. 

The effects of the radiation force and gas pressure on 
the acceleration of winds are summarized in figures 3 and 
4. 

In figure 3 we show the wind final velocity (3 S (thick 
solid curve), the heights, z s and z c , of the wind top and 
the critical point (dashed ones), and the radiation pres- 
sure P c at the critical point (dotted one), as a function 
of the radiative flux F c at the critical point. The other 
parameters are fixed as f = 3, 8 = 1, t c = 1, J = 0.1, and 
c T = 0.3. 

As is seen in figure 3, the wind final velocity increases 
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with F c . The Lorentz factor 7 S of the wind final velocity 
is well fitted by 



7 S ~ 1.07+ 1. 50F C . 



(24) 



Roughly speaking, this is understood as follows. In the 
present treatment under special relativity and pseudo- 
Newtonian potential, there is no rigorous energy integral. 
However, the pseudo total energy along the wind, 



Jc 2 
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(25) 



is approximately conserved within the error of 10%. If we 
denote the gravitational potential by 4>, this pseudo energy 
conservation is written in the non-dimensional form as 
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(26) 



1.0 

0.8 
0.6 
0.4 
0.2 




- 


1 

- 








- -, Zs 











0.1 



0.2 



j_ 



"f(l + 4>) + F 



J_ 

8 2 



7(l+0) + F 



(27) 



Hence, in general 7 S linearly depends on F c . 

In figure 4 we show the wind velocity f3 s (thick solid 
curve), the heights, z s and z c , of the wind top and the 
critical point (dashed ones), and the radiation pressure 
P c at the critical point (dotted one), as a function of the 
isothermal sound speed 6t- The other parameters are 
fixed as f = 3, 6 = 1, r c = 1, J = 0.1, and F c = 0.1. 

As is seen in figure 4, the wind final velocity increases 
with ct- The Lorentz factor 7 S is now fitted by 
,-,2 



l + 0.3c T + 1.364 



(28) 



This is also understood by a rough energy conservation 
discussed above. In this case, 7 S linearly depends on 7 C , 
and 7 C = 1/^/1 — c\, therefore there appears c\ term. 

In addition to transonic solutions, an example of super- 
sonic solutions is shown in figure 5. In figure 5 we show 
the flow velocity j3 (thick solid curve), the radiative flux 
F (dashed one), the radiation pressure P (dotted one), 
and the optical depth r (solid curve), as a function of 
the height z for j = 3, 6 = 1, r = 2, F = 1, /3 Q = 0.6, 
6t = 0.57, and J = 0.1, where the subscript means the 
values at z = 0. 

In this example of supersonic winds, the wind final 
speed is almost the speed of light, since both the initial 
flux and isothermal sound speed are large. 

4. Concluding Remarks 

In this paper we have examined the relativistic radia- 
tion hydrodynamical winds from a luminous accretion disk 
in the relativistic regime of (v/c) 2 . The wind is assumed 
to be steady, vertical, and isothermal. Using a velocity- 
dependent variable Eddington factor, the basic equations 
can be numerically solved without meeting the patholog- 
ical singular point at v = cj\f% to reach the relativistic 
regime. 

For less luminous cases, vertical disk winds are tran- 
sonic types passing through saddle-type critical points, 



Fig. 3. Wind final velocity f3 s at the wind top (thick solid 
curve), heights, z s and z c , of the wind top and the critical 
point (dashed ones), and radiation pressure P c at the critical 
point (dotted one) , as a function of the radiative flux F c at the 
critical point. For heights, z/100 is plotted. The parameters 
are f = 3, 8 = 1, r c = 1, J — 0.1, and c+ = 0.3. 

and the wind final speeds are 0.4-0.8 c for typical pa- 
rameters. For luminous cases, on the other hand, disk 
winds become supersonic types without passing through 
any critical points, and the wind final speeds becomes 
on the order of c. The boundary between the tran- 
sonic and supersonic types is located at around F c ~ 

0.1(e + p)/( P c 2 )/7c. 

In usual standard accretion disks, the local luminosity 
is sub-Eddington. Hence, for proton-electron normal plas- 
mas the above condition is safficicntly fulfilled, and tran- 
sonic winds driven by radiation and gas pressures would 
blow off. For electron-positron pair plasmas, on the other 
hand, the local luminosity becomes super-Eddington, and 
pair winds would supersonically blow off. 

The maximum attainable velocity of the wind can be 
roughly estimated by the energy conservation. In the 
Schwarzschild space-time, without any energy source ex- 
cept for gravity, the total energy of the radiation hydro- 
dynamical flow along the streamline, 



s + p 

Jc 2 — 5-7V500 + gooF = E (const.), 



(29) 



where goo = r g/-Rj is conserved. If all of the thermal 
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Fig. 4. Wind final velocity /3 S at the wind top (thick solid 
curve), heights, z s and z c , of the wind top and the critical 
point (dashed ones), and radiation pressure P c at the criti- 
cal point (dotted one), as a function of the isothermal sound 
speed ct- For heights, 5/100 is plotted. The parameters are 
f = 3, S = 1, r c = 1, J = 0.1, and F c = 0.1. 



energy of the gas and the radiative energy in the initial 
state is converted to the bulk energy of the gas, the final 
velocity would become its upper limit: 



Jc 2 7o 



.9oo F 



(30) 



If, further, the gas is at rest in the initial state and in the 
virial state, [(s + p) / (pc 2 )]y/goo\ z= o = 1, then we finally 
have the maximum attainable velocity as 

ffoo^lo -,,{-, 1 \ F 



7o 
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Jc 2 



1 
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(31) 



0.1, 



For example, in the present typical case of f = 3, F • 
J = 0.1, we have = 1.666, or (3^ = 0.8. 

Hence, from the view point of energetics, the final veloc- 
ity becomes higher and higher, as the ratio Fq/ J becomes 
large; the radiative flux is so high or the mass-loss rate 
is so low. In other words, if the small amount of the gas 
gains the large amount of the radiation energy, the final 
bulk velocity becomes high. In general, however, it may 
be difficult for such a re-distribution of energy to take 
place, and the ratio Fq/ J would be on the order of unity. 

Thus, we conclude that the maximum attainable veloc- 
ity of accretion disk winds emanating from the disk inner 
region would be 0.8 c or so, as long as it consists of normal 
plasmas. This conclusion is consistent with the current 
theoretical works refered in the introduction, those ex- 
plain mildly relativistic jets of 0.26 c (7 = 1.04) in SS 433 
to highly relativistic jets of 0.92 c (7 = 2.55) in several 
microquasars. 

For ultra-relativistic jets of 0.99 c (7 = 10) supposed in 
several active galactic nuclei or extremely relativistic jets 
of 0.9999 c (7= 100) expected to gamma-ray bursts, it may 
be necessary some other processes, including an extra- 
ordinary re-distribution of energy, pair dominant plasmas, 
an energy deposition from other energy sources, such as a 
nuclear one via neutrino, and so on. 




Fig. 5. Flow velocity f3 (thick solid curve), radiative flux F 
(dashed one), radiation pressure P (dotted one), and optical 
depth t (solid curve), as a function of the height z. The 
parameters are f = 3, 5 = 1, to = 2, Fo = 1, /3o = 0.6, ct = 0.57 
and j = 0.1. 
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